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Abstract. We show that a positive signal in a dark matter (DM) direct detection 
experiment can be used to place a lower bound on the DM capture rate in the Sun, 
independent of the DM halo. For a given particle physics model and DM mass we 
obtain a lower bound on the capture rate independent of the local DM density, ve¬ 
locity distribution, galactic escape velocity, as well as the scattering cross section. 
We illustrate this lower bound on the capture rate by assuming that upcoming di¬ 
rect detection experiments will soon obtain a signihcant signal. When comparing the 
lower bound on the capture rate with limits on the high-energy neutrino flux from the 
Sun from neutrino telescopes, we can place upper limits on the branching fraction of 
DM annihilation channels leading to neutrinos. With current data from IceCube and 
Super-Kamiokande non-trivial limits can be obtained for spin-dependent interactions 
and direct annihilations into neutrinos. In some cases also annihilations into rr or 
bb start getting constrained. For spin-independent interactions current constraints are 
weak, but they may become interesting for data from future neutrino telescopes. 
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1 Introduction 

From gravitational effects we know that dark matter (DM) constitutes a significant 
fraction of the energy density of the Universe. Among the most promising ways to 
search for non-gravitational manifestations of DM particles are direct detection (DD) 
experiments, which are looking for the scattering of DM particles from the galactic 
halo in underground detectors [1-8], and neutrino telescopes looking for high-energy 
neutrinos from the annihilations of DM particles in the Sun [9-13]. The latter signal 
emerges from the capture of galactic DM particles in the gravitational potential of the 
Sun after loosing enough energy in a DM-nucleus scattering event in the Sun. Hence, 
the capture rate of DM in the Sun is determined by the scattering cross section of the 
DM particle on nuclei, the same process which provides the signal in DD experiments. 
The resulting neutrino flux from the Sun will also depend on the annihilation channels 
of the DM particle. Therefore, viable information on DM properties can in principle 
be obtained by comparing the two signals. 

However, the DD signal and the DM capture in the Sun depend on different 
parts of the DM velocity distribution. While DD experiments are sensitive to DM 


1 




particles with velocity larger than a certain minimal velocity (which depends on the 
mass of the nuclei in the detector and its energy threshold as well as on the DM 
mass), the DM capture in the Sun is sensitive to values below a certain maximum 
velocity, above which capture of DM particles is kinematically forbidden. Therefore, 
in order to explore the complementarity of the two signals, it is common to adopt 
specihc DM velocity distributions, for instance the so-called Standard Halo Model 
(SHM) consisting of a truncated Maxwellian distribution, see refs. [14-21] for a very 
incomplete list of examples of this approach. However, the properties of the DM 
velocity distribution as well as the local DM density are plagued with large uncertainties 
and halo-independent methods are desirable to draw robust conclusions about possible 
signals. The impact of variations of halo properties on the neutrino signal have been 
studied for instance in refs. [22, 23] and in the context of the neutrino/DD comparison 
in ref. [24] . The authors of ref. [25] investigated the potential to extract DM parameters 
from a combination of data from DD experiments and from a neutrino signal based on 
a polynomial parameterization of the DM velocity distribution, whose parameters are 
htted together with the DM parameters. 

In the present paper we develop a completely halo-independent method to com¬ 
pare a signal from a DD experiment with a neutrino signal from the Sun. We show 
that from a precise measurement of the DD nuclear recoil spectrum a halo-independent 
lower bound on the capture rate in the Sun can be derived. It is based on the overlap 
region in velocity space and therefore does not require any assumptions about the halo 
properties. Our bound extends the halo-independent methods developed in the context 
of DD [26, 27] to the capture rate in the Sun. 

The remainder of this paper is organized as follows: We review the phenomenology 
of DD in section 2 and of the DM induced neutrino signal from the Sun in section 3. We 
then discuss the relation of direct detection to the capture rate in the Sun in section 4, 
where the central result of the paper (the lower bound on the capture rate) is given 
in section 4.2. In section 5, we apply the bound to mock data from future direct 
detection experiments and compare them to the upper bounds from the IceCube and 
Super-Kamiokande neutrino telescopes. We also comment on the importance of nuclear 
form factor uncertainties as well as on the ratio of the neutron and proton couplings. 
We summarize and give our concluding remarks in section 6. In appendix A, we discuss 
how to use an annual modulation signal in a DD experiment to provide the lower bound 
on the capture rate, and apply our results to the DAMA signal. In appendix B we show 
how to use our results in the case of more general scattering cross sections, beyond 
contact interactions. 

2 Dark matter direct detection 

In this section we review the relevant expressions for DD of dark matter [28] . We focus 
on elastic scattering of DM particles x with mass off a nucleus with mass number A 
and mass ruA, depositing the nuclear recoil energy En. The differential rate (measured 
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with being the local DM mass density and Vm is the minimal velocity of the DM 
particle required for a recoil energy Er: 
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where p^A is the reduced mass of the DM-nucleus system. The function fdet{.v,t) 
describes the distribution of DM particle velocities in the detector rest frame, with the 
normalization 
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and we define the angular averaged velocity distribution function / by 


j{v) = Jdttfiv.n), 


(2.4) 


where dO. = d cos 6 dcp. The velocity distributions in the rest frames of the detector, 
the Sun and the galaxy are related by /det(D t) = fsun{v + Ve(t)) = fga,i{v + Vs + Ve(t )), 
where Ve{t) is the velocity vector of the Earth relative to the Sun and Vg is the velocity 
of the Sun relative to the galactic frame. The revolution of the Earth around the Sun 
encoded in Ve{t) leads to an annual modulation of the DD signal [29, 30]. 

To be specihc, in the following we will concentrate on spin-independent (SI) and 
spin-dependent (SD) scattering from a contact interaction. This implies that the dif¬ 
ferential scattering cross section daA{v)/dER scales as 1/n^. Our results can also be 
generalized to other v dependences, as shown in appendix B. For SI contact interactions 
with equal DM couplings to neutrons and protons the cross section becomes 
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where is the total DM-proton scattering cross section at zero momentum transfer, 
p^p is the DM-proton reduced mass, and Ea{Er) is a nuclear form factor. For SD 
interactions a similar formula applies with no enhancement and a different form 
factor. 

Then the event rate eq. (2.1) becomes 


E-iER, t) = A^E\{Er) r]{vm, t ), ( 2 . 6 ) 

^For detectors with several nuclei the total rate is the sum of the rates in all nuclei, i.e., TI{Er, t) = 
'AjaEa(ErP), but usually one nucleus gives the dominant contribution to the rate for a particular 
DM mass. 
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where Vm is considered as a function of according to eq. (2.2) and we have defined 

fj{vm,t) = Cr]{vm,t) with T]{vm,t) = / dvvfdetiv,t) and C = ^ \ . (2.7) 

Jvrr. ^ruxpcp 

Notice that fj depends only on the DM properties and the DM halo, but is independent 
of the characteristics of the experiment. For fixed DM mass, one can translate the event 
rate in Er space into Vm space, and 'fj{vm, t) then has to be the same for any experiment. 
This is the basis of the halo-independent methods developed in refs. [26, 27] and used 
extensively to compare results of different DD experiments, see, e.g., refs. [31-42]. 

For a specihc detector the number of DM induced events in an energy range 
between Ei and E 2 is given by 

poo 

N[Ei,E2] = MT / dEjiE^{ER)G[Ei,E2]iER)fj{vrmt), ( 2 . 8 ) 

Jo 

where M and T are the detector mass and exposure time respectively, and Giei,e 2 ]{Er) 
is the detector response function describing the probability that a DM event with true 
recoil energy Er is reconstructed in the observed energy interval [Ei , E 2 ], including 
energy resolution, energy dependent efficiencies, and possibly also quenching factors. 

In section 4.2 we will assume that a positive signal is observed in a direct detection 
experiment. In this case, the angular averaged distribution f{v) times the constant C 
can be extracted from the data (modulo experimental resolutions and uncertainties). 
Using eqs. (2.6) and (2.7), we hnd [43] 
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where Er is considered as a function of v = Vm according to eq. (2.2), depending on 
the DM mass. 


3 The neutrino signal from DM annihilations in the Sun 


In this section we briefly review the relevant expressions for the neutrino signal from 
the Sun. Due to the scattering of DM particles with the nuclei in the Sun they may 
lose energy and become gravitationally bound to the Sun. Their annihilation products 
can produce neutrinos with energies comparable to the DM mass, detectable at Earth 
[44-46]. 

The capture rate of dark matter particles is given by (see, e.g., refs. [46, 47]) 


Csun — 



dv fiy) V w D^(tc, r), 


(3.1) 


where, for SI interactions, the sum over A goes over all element abundances in the 
Sun up to nickel, while for SD interactions only hydrogen is relevant. We use the 
short-hand notation /sun = / and the angular averaged velocity distribution function 
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is defined in eq. (2.4). The fnnction f{v) refers to the velocity distribntion at inhnity, 
whereas w is the DM velocity inside the gravitational potential of the Snn at the radial 
distance r from the solar center, given by with Mesc(^) being the 

escape velocity from the Sun depending on the location r. The quantity flA{,w,r) is 
the rate with which a DM particle with velocity w will be gravitationally captured by 
scattering on the nucleus A in a spherical shell at radius r with thickness dr\ 


r) = w 

mA 


i(^) 


dEn^iw) 
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where Pa(^) is the mass density of the element A in the Sun {pa/^a being the number 
density). For the numerical calculations we use the standard solar model from ref. [48]. 
Here, Er is the recoil energy of the nucleus after the scattering. For a given DM velocity 
w there is a minimal and maximal nuclear recoil energy, E-cAmiw) and E m»^ (w), such 
that the DM particle gets trapped, i.e., its velocity after the scattering is less than the 
local escape velocity. Writing them in terms of the velocity v at inhnity one has 

Emin = , Emax = -^{v^ + mL(^)) ) (3-3) 

where the former constraint is the requirement for the DM particle to be captured and 
the latter is based on the maximal energy transfer allowed by kinematics. In hgure 1 we 
show Emin and Emn^ir) versus velocity v for scattering on hydrogen for two different 
DM masses, = 10 and 100 GeV, and for the two extreme values of the escape 
velocity in the Sun: 618kms“^ at the surface and 1381 km s“^ in the center of the Sun. 
One can see that F^min and F^max(’") cross, which dehnes the maximum velocity the 
DM particles can have in order to be trapped in a single interaction, n^oss(^)- Froiii 
eq. (3.3) we hnd 


^ Am Amy, 
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Therefore, ncross(’") decreases with DM mass m^, and the capture rates will decrease 
and eventually vanish for large DM masses. For SD interactions, where only hydrogen 
is relevant, we can use mp my. and 


/ ‘fYl 

<ossir) ~ —Uesc{r) or ^ 2mpU^^^{r) . (3.5) 

V 

Typically, for heavier nuclei n^oss(^)) which is relevant for SI interactions, is larger than 
'Across(^)) relevant for SD interactions. 

The cross section daA{w)/dER in eq. (3.2) is the same as the one entering in 
the event rate for DD experiments. Again, restricting to contact interactions^, we use 
eq. (2.5) for SI interactions and an analogous expression for SD. For the nuclear form 

^For other types of interactions see appendix B. 
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Figure 1. For hydrogen (SD), we show in blue the lower bound on the energy for DM 
capture in the Sun, E m\u , and in red the energy upper limits, for the two extreme 

escape velocities, versus the velocity v, for = 10 (100) GeV, left (right). The points were 
they cross (indicated by the vertical dotted lines) are the maximum velocities of the DM 
particles, vlross{r). 


factors of the elements in the Sun we use the approximation F‘^(Er) ~ e with 

Ea = 3/{2mARA) and Ra = [0.91(m^/GeV)^/^ + 0.3] fm [46]. Hence, we obtain 

_ P^Sun r^cross _ r^max(^) 

Gsun = 47r C ^ drr'^pAir) dvf{v)v / Fl{ER)dER (3.6) 

4-^0 Jo '^^min('h’) 


for the capture rate in SI interactions. The coefficient C is defined in eq. (2.7) and 
contains the DM-nucleus scattering cross section. For SD interactions only hydrogen 
is relevant, i.e., the sum contains only one term with H = 1, and the form factor is 
trivial, F^{Er) = 1. 

If equilibrium between DM capture and annihilation in the Sun is reached, the 
hnal annihilation rate is independent of the annihilation cross section and is given by: 

rSun=2^Sun- (3.7) 


For equilibrum to occur, the equilibration time teq must be smaller than the age of the 
Sun, i.e., Gq < tsun 4.5 Gyr, where (see for instance ref. [19]) 
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\/Gsun ^Sun 
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(3.8) 

Here Hsun is the annihilation rate in the Sun and {av) the thermal average of the 
annihilation cross section. As we will see in section 5 for mock data and in appendix A 
for DAMA, our values of the lower bound on the capture are safely above > 10^^ s”^ in 
all cases except for SI interactions in xenon, due to the small scattering cross section 
assumed for the mock data, cxsi = lO-'^^ QjqqS g-i_ equilibrium may not be 

reached for annihilation cross sections smaller or equal than the freeze-out one. Also 
we note that in case of p-wave annihilations the annihilation cross section today can 
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be much smaller than the thermal freeze-out value assumed in eq. (3.8) and depending 
on the capture rate equilibrium may or may not be reached. 

Notice that we are neglecting evaporation, which is justihed for > 3 GeV 
[45, 46, 49]. The neutrino flux at the Earth from the annihilation channel / with 
branching ratio BRj is 


M = rr rsu„ dNl 
dEy ^ An dEy ’ 


(3.9) 


where d is the Sun-Earth distance and dNl/^Ey is the neutrino spectrum per annihi¬ 
lation of flavour / that reaches a distance of 1 AU, which needs to take into account 
flavour transitions (neutrino oscillations, MSW effect), absorption and regeneration, 
see, e.g., refs. [50, 51]. 


4 Relating DM direct detection to the DM capture rate in 
the Sun 

We now want to relate a DD signal to the DM capture rate in the Sun in a halo- 
independent way. We assume that the DM velocity distributions relevant in the two 
cases are the same. This implies two important consequences: 

1. We neglect the small velocity of the Earth with respect to the Sun and adopt the 
approximation for / (dehned in eq. (2.4)) 

/det(^^) ~ hnn{v) = f (v) . 

Hence we use that Vg ~ 29km/s -C Vm for typical values of Vm- This also implies 
that we are ignoring the (small) annual modulation signal in the direct detection 
rate (we will comment on this in appendix A). 

2. DM direct detection samples the DM distribution today, while for the DM capture 
in the Sun the velocity distribution on time scales relevant for the equilibration 
of capture and annihilations is relevant. We will assume that the DM velocity 
distribution is constant on those time scales and that the same f{v) applies for 
the neutrino capture and direct detection. Similarly, we also assume that the 
energy density is constant on times scales relevant for equilibration and equal to 
the current one. 

4.1 The overlap in Vm 

From inspection of the relations given in sec. 2 we see that direct detection is sensitive 
to high DM velocities. Let us denote the threshold energy of a given experiment 
by Ethr- If we ignore the hnite energy resolution of the experiment, this dehnes a 
threshold velocity via eq. (2.2): Uthr = 'ym(T^thr)-^ Hence, this experiment is sensitive 

^In the case of a finite energy resolution the actual threshold velocity would be somewhat smaller 
than the value corresponding to the nominal threshold energy of the experiment due to the recon¬ 
struction of events below the threshold at higher energies. 
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Figure 2. In blue we show the minimum velocity Vthr(-£'thr) probed in a direct detection 
experiment versus m^, assuming different threshold energies Fthr and using a Ge (left) and a 
Xe (right) target. In red we show the maximum velocity relevant for DM capture in the Sun 
for scattering on hydrogen, Ucross, for the two extreme values of Uesc = 1381 kms“^ in the 
centre of the Sun and Vesc = 618kms“^ at the surface. The shaded area shows the overlap 
region assuming scattering in the centre of the Sun. The horizontal black line indicates 
approximately the galactic escape velocity in the detector rest frame. 


to DM velocities v > Uthr- In a realistic halo there will be a maximal velocity of DM 
particles, set by the escape velocity from the galactic gravitational potential, Uesc- Since 
the precise value of this escape velocity is uncertain (see, e.g., ref. [52] and references 
therein) we want to keep the discussion as independent of this argument as possible. 

In contrast to direct detection, eq. (3.6) shows that only DM velocities v < 
contribute to the capture rate in the Sun. Hence, direct detection is sensitive to DM 
velocities above a certain value given by Uthr, whereas for DM to be captured in the 
Sun velocities below the value set by are relevant. In order to relate the two 

phenomena we need to consider the overlap regions in velocity space. If v^oss < 
there is no overlap and direct detection and DM capture decouple. In this case no 
statement can be made halo-independently and a connection can be established only 
by referring to some a-priori assumptions about the DM velocity distribution. 

In hgure 2 we show with blue curves the minimum velocity probed in direct 
detection experiments versus the DM mass, using Ge (left) and Xe (right) for different 
examples of threshold energies. Those curves can be compared with the maximal 
velocity relevant for the DM capture in the Sun, Ucross, which is shown with red curves 
for scattering on hydrogen for the two extreme values of the solar escape velocity Uesc 
corresponding to the values in the centre and at the surface of the Sun. From the 
hgure we observe that for a large range of DM masses there is an overlap of the DM 
velocities probed by direct detection and solar capture, i.e., Vthr < Vcmss, at least for 
scattering events which occur near the centre of the Sun. A similar hgure can be 
found in ref. [25], where it is also shown that Across is even larger for heavier nuclei 
relevant for SI interactions and therefore the overlap is also enhanced with respect to 
SD interactions. 

In the following we are going to make use of the overlap region in velocity space in 










order to compare results from direct detection experiments to neutrino searches from 
the Sun without specifying any halo properties apart from the two assumptions stated 
at the beginning of this section. In particular, in section 4.2 below we are going to 
assume a positive signal in direct detection and derive a halo-independent lower bound 
on the capture rate in the Sun. We will use the fact that from a precise measurement of 
the nuclear recoil spectrum the angular averaged distribution f{v) times the constant 
C can be extracted from the data via eq. (2.9), and we can make use of the signal 
precisely in the overlap region in velocity space. 

Before we proceed with this, let us briefly comment on the case of a positive signal 
from neutrinos. In this case no direct information can be obtained about f{v), and the 
signal can come from any DM velocity below Ucross- For a positive neutrino signal no 
general lower bound for the direct detection rate can be established halo-independently, 
since the capture may happen entirely from velocities below the threshold of the DD 
experiment. If no signal is seen in direct detection at the relevant level, such a situation 
could point to a halo dominated by low velocity DM particles (in the solar frame), for 
instance a co-rotating dark disk [22]. Alternatively this may indicate more exotic 
particle physics such as self-interacting DM, which enhances the capture rate by DM- 
DM scattering, independently of the DM-nucleus scattering cross section [53]. Another 
example of a modihed relation between the DD signal and capture in the Sun is inelastic 
DM scattering [54-56]. These cases will be studied in a future work. 

4.2 A lower bound on the capture rate from a positive direct detection 
signal 

We can derive a lower bound on the capture rate, eq. (3.6), by using that for any positive 

A A 

function H{v) > 0 we have that if(n) dn > H{v)dv for Ucross > 

Therefore we have 

_ P^Sun ^"^cross _ 

Csun > dvr ^ / drr^pA{r) dv f{v) v J^Aiv^r), (4.1) 

^ 0 j fJthr 


where we dehned 


J^Aiy.r) 


rEuia^iv) 


Fl{ER)dER, 


(4.2) 


with Erainiv) and F^niax(t') given by eq. (3.3). Notice that EA{v,r) and 

eq. (3.4), depend on r via Uescij")- 

Let US assume now that a signihcant direct detection signal has been observed. 
We can then use eq. (2.9) to obtain f{y) from the data and insert it in eq. (4.1): 


Csun>47r5^A2 
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= 47r 


"Rsu 


'•^Su 


drr‘^PA{r) 


drr‘^pA{r) 


‘^cross 


dv 


'Vthi 
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EA{v,r) 


vivthr)J^A{vthr,r) + / dv p{v) E'a{v, r) 


'd^thr 


(4.3) 
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where in the last line we integrated by parts, with = 0 = 

dJ^A{v,r)/dv. For SD interactions only hydrogen is relevant and we have Fii{Eji) — 1, 
and 




(4.4) 


The direct detection signal has to be precise enough for fi{v) to be extracted from the 
observed energy spectrum via eq. (2.9) with sufficient precision. Either the derivative 
or the function fj{v), including its value at the experimental threshold, have to be 
determined, including unfolding of experimental resolutions and backgrounds. This 
will require a signihcant number of events such that an accurate spectral analysis can 
be performed. If those conditions are met, eq. (4.3) provides a lower bound on the 
capture rate in the Sun without specifying the DM velocity distribution, the galactic 
escape velocity, the scattering cross section or the local DM density. This lower bound 
is the central result of this paper and we will illustrate it numerically in the following 
section for possible future signals in DD experiments. 

5 Numerical examples 

5.1 Mock data for direct detection 

The halo-independent bound on the DM capture rate in the Sun derived in the previous 
section will be useful once a clear and highly significant signal in a direct detection 
experiment has been observed.^ In the following we will assume that the DM-nucleon 
scattering cross section is just below the current limits [5, 8], which will allow upcoming 
experiments to obtain a signihcant signal. As representative examples we will consider 
a future xenon based experiment [57-59] as well as a germanium detector [60, 61]. 
For calculating mock data we adopt the conventional Maxwellian velocity distribution 
(SHM) with V = 220 km/s, truncated at the escape velocity of Vesc = 544 km/s, and 
we assume a local DM density = 0.3 GeV/cm^. The velocity of the Sun in galactic 
coordinates is (10, 233, 7) km/s. 

For the xenon experiment we adopt a threshold of 5 keV (we also show the effect 
of a reduced threshold of 3 keV below) and we take natural abundances of the isotopes 
with spin, ^^®Xe (26.4%) and ^^^Xe (21.2%). We use the values of 10“^^ cm^ for the 
SI cross section and 2 ■ 10“'^° cm^ for the SD cross section, with equal couplings to 
protons and neutrons in both cases. Assuming = 100 GeV, for an exposure of 

'*’The CDMS collaboration reports 3 candidate events from their data with a silicon target, with 
a 0.19% probability for the known-background-only hypothesis when tested against the alternative 
DM-l-background hypothesis [7]. In order to apply our bounds a detailed spectral measurement 
is required and not enough information can be extracted from the 3 observed events in CDMS-Si. 
However, motivated by this potential signal, we take DM parameter values similar to those preferred 
by CDMS-Si events for our example of mock data in germanium. The annual modulation signal 
reported by the DAMA/LIBRA collaboration [1] will be discuss in appendix A. We emphasize that 
both signals are either excluded or in strong tension with several other experiments [2-6, 8] (see for 
instance refs. [35, 42[ for halo-independent analyses). 
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(GeV) 


Figure 3. Lower bounds on the DM capture rate in the Sun for xenon (left) and germanium 
(right) DD mock data as described in the text, both for SI (red curves) and SD (blue curves) 
interactions. For illustrative purposes, in the Xe SD case we show the effect of changing 
the default threshold energy of Ethr = 5 keV (solid) to 3 keV (dashed). The dotted curves 
correspond to the actual capture rate. To calculate the DD mock data as well as the capture 
rate we assume the standard halo model. 


1 ton yr at 100% efficiency abont 154 (267) events would be observed in the energy 
range 5 — 45 keV for SI (SD) case. 

We also generate mock data for a future germanium experiment, with a threshold 
of 1 keV, focusing on low DM masses. Assuming a DM mass = 6 GeV and cross 
sections of cxsi = 5 ■ lO”"^^ cm^ and ctsd = 2 ■ 10“"^° cm^ (with equal couplings to protons 
and neutrons) we would obtain about 1.5 x 10^ (2-3) events for SI (SD) interactions in 
the energy range 1-10 keV for an exposure of 100 kg yr with energy resolution of 30%. 
While the SI case will allow for a high statistics reconstruction of the event spectrum, 
this will not be possible for the SD case with the assumed parameters. This follows 
from the missing enhancement with respect to SI interactions, as well as the small 
abundance of only 7% of the relevant Ge isotope with spin. Nevertheless we are going 
to include SD interactions in the following discussion, keeping always in mind that this 
will require larger exposures or a different target nucleus with better sensitivity to SD 
interactions. In addition, the bounds are stronger for larger DM masses, but we want 
to provide an illustrative example for the small mass region. 

In the following analysis we neglect the energy resolution, effects of binning the 
data, possible contamination with background, experimental errors, and nuclear form 
factor uncertainties (we comment on form factors later in subsection 5.4). Hence, we 
assume that fjiym) can be extracted from the observed nuclear recoil spectrum (for a 
specihed m^). This idealized analysis suffices to illustrate the power of our bounds. 
Once applied to real data an appropriate statistical analysis will have to be performed. 

In hgure 3 we show the lower bound on the capture rate for the mock data of future 
DD xenon (left) and germanium (right) experiments as described above, both for SD 
and SI interactions and compare them to the true capture rate assuming the SHM. In 
the case of Ge we focus on the low DM mass region, i.e., < 30 GeV, keeping in mind 

that the bounds are stronger (i.e., closer to the true capture rate) for > 20 GeV. 
The bounds as well as the capture rate are shown for the “true” DM mass, which has 
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been used to calculate the mock data. We see that our bounds are strong in a large 
portion of parameter space for the xenon experiment: 20 < < 1000 GeV for SD, 

and for > 50 GeV for SI. 

5.2 Limits from neutrino telescopes 

In order to compare the lower limit on the DM capture rate in the Sun from a DM direct 
detection signal with the upper bounds on the neutrino rate from neutrino telescopes 
we proceed as follows. We assume equilibrium between capture and annihilations and 
the neutrino-induced muon rate in neutrino telescopes will then depend on the specihc 
annihilation channel, which depends on the particle physics model for DM. In partic¬ 
ular cases, the high-energy neutrino rate may be strongly suppressed by the available 
annihilation channels, for instance for annihilations into or u,d,s quarks (see, 

however, ref. [62] for higher order effects). Hence, translating a given upper bound on 
the neutrino-induced muon rate into a bound on the capture rate depends on the DM 
annihilation channel, see section 3. 

Below we are going to compare the lower bounds on Gsun to the limits from Ice- 
Gube (IG) and Super-Kamiokande (SK). For IG we use the upper bounds of ref. [11], 
where results are given directly as upper limit on the capture rate for various annihila¬ 
tion channels as a function of the DM mass. We will show results for two cases, namely 
annihilations into bb + rr and into WW + rr (we keep particle/antiparticle notation 
implicit for annihilation products). For SK there are no limits on the capture rate 
available. Therefore, for annihilation channels into bb and rr in the low DM mass re¬ 
gion (up to 200 GeV) we extract the upper limits on the capture from the upper limits 
on the scattering cross section given by the most recent SK results [13]. For the large 
mass region in these channels, and also for direct annihilation into neutrinos (r'^^r'^), 
we use the limits on the capture rate calculated in ref. [63] (Tab. II) based on SK data 
from ref. [10]. The quoted bounds at the 90% GL from IG and SK apply assuming 
that annihilations proceed with 100 % branching ratio into the indicated channels. 

5.3 Comparison of direct detection and neutrino data 

In order to apply the lower bound on the capture rate derived in section 4.2 one has 
to specify the DM mass and the couplings to neutrons and protons. For the moment 
we will restrict the analysis to equal couplings to neutrons and protons, and we will 
come back to this point in section 5.4. Regarding the DM mass, in general it cannot be 
extracted from a DD signal without referring to a specihc DM halo model. Therefore, 
we do the analysis without assuming that the DM mass is known. We calculate mock 
data for a hxed “true” DM mass, but then we apply the lower bound on Csun as a 
function of (different from the “true” value). This procedure resembles the situation 
we would face in case of applying the method to real data. If the mass was known 
from some other data (e.g., an observation at LHG or a 7 line signal from indirect DM 
searches) one would of course perform the analysis only for that DM mass. 

In hgure 4 we show the lower bounds on the capture, assuming a true (but un¬ 
known) dark matter mass of = 30 GeV (upper panels), 100 GeV (middle panels), 
and 300 GeV (lower panels) for a xenon DD experiment and compare them to the 
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Figure 4. Left: Lower bounds on the capture rate from a future xenon DD experiment 
for SI (red) and SD (blue) interactions compared to the 90% CL upper limits for the com¬ 
bined annihilation channels WW^tt and bb,TT from IceCube (IC, green curves) [11] and the 
channels rr, bb from Super-Kamiokande (SK, black curves) [10, 13, 63]. Right: Upper 
bounds on the branching ratios versus dark matter mass for SD interactions. To calculate the 
DD mock data we assume a “true value” for the DM mass of niy. = 30, 100, 300 GeV in the 
top, middle, bottom panels, respectively, and values of 10“^^ cm^ for the SI cross section and 
2 • 10“^^ cm^ for the SD cross section, with equal couplings to protons and neutrons in both 
cases. Assumptions about the mock data for the DD experiment are given in section 5.1. 

limits from IC and SK. While for SI interactions the lower bonnd from DD would be 
consistent with the limits from IC and SK, we see that for SD interactions tension 
arises if DM annihilates into neutrinos or into tt/WW. In the right panels we show 
the ratio of the upper limit on Csun to the lower bound. This can be interpreted as an 
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Figure 5. Left: Lower bounds on the capture rate from a future germanium DD experiment 
for SI (red) and SD (blue) interactions compared to the 90% CL upper limits for the combined 
annihilation channels WW^tt from IceCube (IC, green curves) [11] and the channels 
rr, bb from Super-Kamiokande (SK, black curves) [10, 13, 63]. Right: Upper bounds on 
the branching ratios versus dark matter mass. To calculate the DD mock data we assume a 
“true value” for the DM mass of m-^ = 6 GeV and cross sections of cisi = 5 • 10“^^ cm^ and 
Usd = 2 • 10 “^° cm^ (equal couplings to protons and neutrons). Assumptions about the mock 
data for the DD experiment are given in section 5.1. 


upper bound on the branching ratio of the corresponding annihilation channel. This is 
conservative, since it assumes that there are no neutrinos from any other annihilation 
channels, which typically will not be the case. 

We see that direct annihilations into neutrinos would be constrained to branching 
ratios at the few % level, with some dependence on the DM mass. Annihilations 
into rr, WW could be constrained at the 10% level. However, note that for true 
< 100 GeV the bound is stronger when applied at a “wrong” DM mass, i.e., for 
values of which are larger than the actual true value. When using the correct value 
of the limits for the assumptions adopted here would be rather weak. 

In hgure 5 we show the lower bounds on the capture (left) and the upper bounds 
on the branching ratios for different channels (right), assuming the low threshold ger¬ 
manium experiment. Since this conhguration is most sensitive at low DM masses we 
take here a “true” DM mass of 6 GeV. In this case SD and SI interactions lead to sim¬ 
ilar lower bounds which potentially can constrain annihilations into channels leading 
to neutrinos. Again we note the feature that bounds get stronger when applied for the 
“wrong” DM mass, in which annihilations into rr become also constraining. 

We note that our comparison of DD and neutrino data is not fully consistent, 
in the sense that we are comparing a possible signal in a future DD experiment with 
current limits from neutrino telescopes. When the potential signals from DD will be 
available, limits from IG and/or SK may have improved and upgraded and/or new 
neutrino telescopes may be in place [64-66]. In this sense our results are conservative, 
since the comparison may become more stringent than what is shown here. 

Let us comment briefly on the annual modulation of the DD signal. We expect 
that using information on the modulation in addition to the unmodulated rate does 
not provide signihcantly stronger bounds on the capture rate. If a future experiment 
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only measures the modulation amplitude but cannot distinguish the unmodulated DM 
rate from background (similar to the DAMA experiment [1]) one can use the results 
of refs. [34, 35] to derive a lower bound on Csum based on the modulation amplitude. 
Those bounds require some additional (modest) assumptions on the halo and they can 
be found in appendix A, where we also apply those bounds for the DAMA signal [1], 
which however is incompatible with other DD limits halo-independently [35]. 


5.4 Uncertainties due to form factors and unknown couplings to protons 
and neutrons 

In this section we discuss the uncertainties associated to nuclear form factors, as well as 
the case of arbitrary couplings to neutrons and protons. Different calculations of form 
factors sometimes lead to large differences, especially for SD interactions and heavy 
nuclei such as xenon, see for instance refs. [67-71] for discussions. In our numerical 
calculations based on mock data no form factor appears, since we assume that the 
“correct” one is used, when extracting the velocity distribution from the observed rate 
via eq. (2.9). If a “wrong” form factor was used in eq. (2.9) it would modify the 
extracted velocity distribution, fextriv), in the following way: 


Cfey,tv{v) = Cf{v) 


KneiEn) 


F2 (Er) 

wrong \ 


vjv) d f \ 

V dv \E^^^^^{Er)) ■ 


(5.1) 


Let us focus on the case of SD interactions (where form factor uncertainties can be 
large) and consider also general couplings to neutrons (a„) and protons (op). In this 
case the form factor can be expressed in terms of the structure functions Sij{ER) as 

F|d(Ur) = (l + «)' S'oo(Uk) + (l - Soi{Er) + (l - nySniER) , (5.2) 


where we dehne n = an/dp and we have absorbed dp in the cross section Now we 
can use this relation together with eq. (5.1) to test both the impact of using a “wrong” 
form factor as well as “wrong” values of k when deriving the bound on the capture 
rate from mock data. 

In order to illustrate those effects we generate mock data assuming SD interactions 
for a neutron-dominated experiment (xenon) and for a proton-dominated one (fluorine). 
We use the same parameters for the mock data as in section 5.1 for both experiments, 
with (TsD = 2 ■ 10“^° cm^, a DM mass of = 100 GeV, and Ethv = 3 keV. For these 
parameter values, we predict 78 events in the energy range [3,10] keV for a fluorine 
experiment with total exposure of 100 kg ■ y. To generate mock data we assume equal 
couplings to protons and neutrons (k = 1 ). 

In hgure 6 we show the lower bound on the capture rate as a function of k, 
illustrating the case of using a “wrong” ratio of neutron to proton couplings when 
analysing the data. For k ~ 1 both F and Xe give similar lower bounds, with the 
difference coming from their minimum velocities at the threshold (uthr) due to their 
different masses. For the xenon experiment, for |fi;| S> 1, the capture rate vanishes, 
since in xenon spin is mostly carried by the neutron. Hence, for IkI ^1 the DD 
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Figure 6. Lower bound on the capture rate for SD interactions from mock data from 
experiments using xenon (red) and fluorine (blue). Mock data is generated for a DM mass 
of = 100 GeV, <7^ = 2 • 10 ““^° cm^, and equal couplings to protons and neutrons {k = 
anlo,p = 1). The lower bound on Csun is shown as a function of |k| (i.e., using a “wrong” 
value of K when using the mock data to calculate the bound), where the solid (dotted) curves 
correspond to k > 0 (k < 0). Dark (light) colours correspond to the form factors from 
ref. [67] (ref. [70]). 


signal is dominated completely by nentrons and since the captnre in the Snn is set by 
protons, Csun becomes snppressed by in that case. On the other hand, the spin 
of the flnorine nnclens is mainly provided by a proton. Hence, for |fi:| < 1 the bonnd 
on the captnre rate becomes independent of the ratio of the conplings, since both DD 
and Osun are controlled by interactions with protons. 

We also illustrate the effect of using different form factors in hgure 6 , comparing 
form factor calculations from ref. [67] (darker colours) and ref. [70] (lighter colours). We 
see that those two examples give similar results in the regime where the interaction is 
“large”, i.e., for neutron (proton) dominated interactions for xenon (fluorine). However, 
in the opposite case those two form factor calculations lead to very different results. 
For |fi:| < 1 (proton domination), a xenon experiment would give very strong bounds on 
the capture rate, since interactions are suppressed for DD, and explaining the assumed 
signal would require a very large scattering cross section. However, in this regime the 
bounds differ by about one order of magnitude between the two form factor models. 
Similarly, for |k| > 1 (neutron domination), the bound on the capture rate from a 
fluorine experiment changes by more than a factor 100 between the two form factor 
calculations. 

Let us stress an important point related to using “wrong” values of k. This dis¬ 
cussion is mostly relevant if data is available only from one experiment (or experiments 
with the same spin structure). If a signihcant signal from xenon as well as fluorine 
is observed (such as assumed in our mock data), then the neutron to proton ratio k 
is essentially determined by the relative strength of the two signals. This emphasizes 
the need of data from complementary targets. A similar discussion will also apply for 
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SI interactions with arbitrary couplings to neutrons and protons. We leave a detailed 
study of SI interactions with general isospin structure for future work. 

In hgure 6 we have changed the form factors both for generating mock data as well 
as calculating the bound on the capture rate, i.e., we have always used the “correct” 
form factor. We have also tested the impact of adopting a “wrong” form factor by 
using eq. (5.1). By comparing the form factors from refs. [67] and [70] in this way 
we hnd in some cases even unphysical negative values for fextr{v)- This means that 
the wrong form factor leads to inconsistent results and one would be able to see from 
the observed spectrum that data are not consistent with that particular form factor 
choice. Hence, we note that form factor uncertainties are important, but once a precise 
spectrum from DD is available (such as necessary for our method to work) we will have 
an additional tool at hand to test different form factor models. 

6 Discussion and conclusions 

We have established a halo-independent framework to relate a signal in a DM direct 
detection experiment to the neutrino rate in neutrino telescopes from DM annihila¬ 
tions in the Sun. Assuming that the DM velocity distribution and the DM density 
are constant on time scales relevant for equilibration in the Sun, we have derived a 
lower bound for the DM capture rate in the Sun in terms of a positive signal in a 
direct detection experiment, see eq. (4.3). If DM capture and annihilation in the Sun 
are in equilibrium we obtain (conservative) upper bounds on branching fractions for 
annihilations in channels involving neutrinos from the comparison of the lower bound 
on the capture rate from a direct detection signal with the upper limits from neutrino 
telescopes. 

The lower bounds are based on the part in DM velocity space which contributes 
to both the capture in the Sun as well as the scattering in direct detection experiments. 
We hnd that for typical threshold energies of direct detection experiments a signihcant 
overlap region exists to apply our bounds. Hence, the lower bounds are independent of 
the velocity distribution, the escape velocity, the cross section or the local DM density, 
although we implicitly assume that these values are such that a direct detection signal 
can be measured. With some additional modest assumptions on the halo properties the 
bound can also be used for an annual modulation signal, as explained in appendix A. 

To illustrate the power of the bounds, we have applied them to mock data from 
future xenon and germanium experiments, assuming that the true DM cross section 
is not too far from the current upper bounds. In such a case the halo-independent 
comparison to present limits from Super-Kamiokande and IceCube leads to non-trivial 
bounds on the branching fraction for direct annihilations into neutrinos, if the scat¬ 
tering is spin-dependent. In some cases also annihilations into rr start getting con¬ 
strained. For spin-independent interactions current constraints are weak. In general 
we note that bounds are stronger for DM masses in the range between 100 to 500 GeV. 
For lower DM masses the lower bound from direct detection as well as the limits from 
neutrino telescopes become weaker. Those results can be found in hgures 4 and 5. We 
note that this method is expected to become more powerful when data from future 
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neutrino telescopes such as IceCube-Pingu, Hyper-Kamiokande, or KM3NET become 
available. 

The halo-independent comparison has to be done for a specihc model for the 
DM-nucleon interaction (and as a function of the DM mass). In our work we mostly 
assumed either SI or SD elastic scattering with equal couplings to neutrons and pro¬ 
tons. Generalization to other types of couplings is straightforward. As an example we 
discussed in section 5.4 the case of SD interactions with arbitrary couplings to neutrons 
and protons. This leads to interesting situations, since only scattering on free protons 
is relevant for the capture in the Sun, whereas for direct detection the sensitivity is 
governed by the spin-composition of the target nucleus. One may imagine a situation 
where the spin of the nucleus is dominated by neutrons, in which case the neutrino 
and direct detection signals largely decouple. In order to exclude such a case it will be 
essential that data from target nuclei with spin carried by protons is available. Another 
case in which a DD signal does not always imply a neutrino signal from the Sun is 
asymmetric DM (see ref. [72] for a recent review). While our bounds on the capture 
rate would still apply, in those models annihilations are often suppressed due to the 
lack of anti-DM, and therefore a possible signal from neutrinos becomes very model 
dependent. 

In our work we always assumed that the differential scattering cross section 
da/dEji is proportional to 1/u^, which is true for DM-nucleon contact interactions, 
but in general this may not be the case for more exotic types of interactions (see for 
instance ref. [37] for generalizing halo-independent methods for direct detection to such 
cases). In appendix B we show that it is straightforward to generalize our lower bound 
on the capture rate to all cases where the velocity and nuclear recoil energy dependence 
of the differential cross section factorizes as da/dEn = g{y)h{Ept). Considering more 
complicated models is beyond the scope of this paper and is left for future work. An¬ 
other particle physics variation for which our results do not apply directly is inelastic 
scattering of the type y -f A^ —)■ y* -f- A^, where the mass difference between y and y* is 
of the order of the DM kinetic energy. This will change the kinematics of the scattering 
and generically for > niy. the capture rate is increased [54-56], see ref. [38] for 
halo-independent considerations in the context of DD. We leave also the generalization 
of the lower bound on the capture rate to this case for future work. 

To conclude, we would like to emphasize that in the presence of a DD signal and 
a signal of neutrinos from the Sun, before proceeding to extract the parameters by 
doing a £t to both, one should hrst check if the lower bounds on the capture derived 
here are fulhlled for some combination of annihilation channels and branching ratios. 
This would provide a halo-independent consistency check for both signals. 
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Scientihc Program “News in Neutrino Physics”, where this work was initiated. 
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A Lower bound on the capture for annual modulation signals 


The change of the velocity of the detector relative to the DM halo due to the Earth’s 
rotation around the Sun leads to an annual modulation of the event rate in a DD 
experiment [29, 30]. We denote the amplitude of the modulation expressed in Vm space 
by Ar^ivm)- Information from a modulation signal can be combined with the halo- 
independent lower bound on the capture rate, eq. (4.3), by using halo-independent 
upper bounds on the annual modulation in terms of the average rate derived in 

ref. [34] and applied in refs. [35, 38]. Those bounds are based on an expansion in the 
Earth velocity Vf,, using that Ve/v is small for v >Vm for typical values of Vm relevant 
for experiments. In such an expansion the time independent rate appears at zeroth 
order, whereas the annual modulation amplitude is of linear order in v^- 

If we assume that the DM velocity distribution is constant on time scales of years 
and constant in space on scales of the Sun-Earth distance, one can derive the bound 


[34] 


ArjiyYo) ^ Ug 


dr] ri{vm) 

dv^a '^m 



(A.l) 


Note however, that the lower bounds on the capture rate in eq. (4.3) require a lower 
bound either on {—dfj/dvm) or on which cannot be obtained in terms of A^(nm) 

alone from eq. (A.l). In principle eq. (A.l) can be re-written as a lower bound on 
{—dfj/dvm) involving both, A^(nm) and Hence, this would require an experiment 

able to determine the modulation amplitude as well as the unmodulated rate. Moreover 
we do not expect a signihcant stronger lower bound on Csun from such a procedure 
for modulation amplitudes which fulfill eq. (A.l) (as they must to be consistent with 
a DM signal). 

Adopting some modest additional assumptions on the halo another bound on the 
modulation amplitude can be derived. If there is only one preferred direction in the 
DM velocity distribution one finds (see ref. [34] for details): 


Arjiym) 


< 


-Ve sin tthalo 



(A.2) 


where a is the angle between the preferred DM direction and the direction perpendic¬ 
ular to the Earth’s orbit. Assuming an observed signal for A^(nm), this provides ns 
with a lower bound on {—dfj/dvm), which can be directly plugged in eq. (4.3) and we 
get: 
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where in the last line we used Ar^{v) = CA^(n). In the following we will adopt the 
most conservative option sincthaio — 1, but notice that in cases where the preferred 
DM direction is set by the direction of the Sun relative to the DM halo (e.g., for a halo 
with a dark disk) one has sinohaio — 0.5 and the bounds scale accordingly. 

Notice that observing annual modulation with typical exposures (and allowed 
cross section values) is extremely hard (see for instance ref. [73]). The DAMA/LIBRA 
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Figure 7. In blue we show the minimum velocity Vthr(-£'thr) probed in a direct detection 
experiment versus m^, assuming different threshold energies Fthr and using a Na (left) and 
a I (right) target. In red we show the maximum velocity relevant for DM capture in the 
Sun, Ucross for scattering on hydrogen for the two extreme values of Uesc = 1381 kms“^ in the 
centre of the Sun and Uesc = 618kms“^ at the surface. The shaded area shows the overlap 
region assuming scattering in the centre of the Sun. The horizontal black line indicates 
approximately the galactic escape velocity in the detector rest frame. 


experiment reports an annual modulation of the signal in their Nal scintillator detector, 
with a period of one year and a maximum around June 2nd with very high statistical 
signihcance [1]. This modulation is strongly disfavoured by other experiments halo- 
independently, both for elastic SI and SD interactions [35] and for inelastic scattering 
[38]. Despite these problems of the DM interpretation of the DAMA modulation signal 
we will use it in the following to illustrate the lower bounds on the capture rate eq. (A.3) 
based on a modulation signal. 

A.l The DAMA modulation signal 

In hgure 7 we show the overlap in Vm space between DAMA and the DM capture 
in the Sun assuming scattering either on Na or on I. Na dominates for DM masses 

~ 20 GeV, whereas iodine is relevant for larger DM masses. For sodium we 
observe a large overlap region. Note that the DAMA threshold of 2 keVee corresponds 
to a recoil threshold of 6.7 keV (Na) and 22 keV (I) for usual quenching factors. This 
implies a relatively small overlap region in the case of iodine, especially for the SD case 
shown in hgure 7. 

To obtain a lower bound on the capture rate from DAMA data we consider a 
binned version of eq. (A.3) and use the observed values of Ari{vm) corresponding to 
the energy bins reported by DAMA, see refs. [34, 35, 38] for details on this procedure. 
We assume either elastic SI or SD scattering on Na or on I with equal couplings to 
protons and neutrons. The resulting lower bounds on the capture rate are shown in 
hgure 8 both for SI and SD, together with the 90% CL upper limits on the capture 
from neutrino telescopes. 

From the left panel we see that for scattering on sodium there is tension between 
the lower and upper bounds implying that for SI DM annihilation into neutrinos and 
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Figure 8. Lower bound on the capture rate from DAMA data assuming scattering on sodium 
(left) or iodine (right) for SI (red) and SD (blue) interactions, using the bound from eq. (A.3). 
Also shown are the 90% CL upper bounds from IceCube [11] for annihilations into WW + tt 
and WW + bb (green curves), and from Super-Kamiokande into bb, rr, and [10, 13, 63] 
(black curves). 


TT (bb) are strongly constrained for > 5 GeV (15 GeV), while SD is excluded for 
all channels. 

For scattering on iodine shown in the right panel of hgure 8 we note the strong 
dependence on the DM mass for SD scattering. This can be understood from hgure 7, 
which shows that only for a small range of DM masses there is overlap in Vm space. 
For SI interactions Ucross is larger (not shown in the plot), leading to larger overlap 
in Vm space, and we observe strong bounds on the annihilations shown in hgure 8 for 
my_ > 35 GeV (and even for > 10 GeV for the neutrino channel). 


B Beyond contact interactions 


In the main text of the paper we have always assumed contact interactions between DM 
and the nucleus, which leads to a 1/n^ dependence of the diherential scattering cross 
section. In more exotic models also other dependences on v and/or Er are possible, 
see for instance refs. [21, 37, 63] in the context of DD and the neutrino signal. Our 
bounds can be generalized in a straightforward way if the dependence of v and 
factorizes. Let us assume that the diherential cross section on the nucleus with mass 
number A can be written in the form 

^=gA(v)hA(EK). (B.l) 

In the conventional case considered in the main text we have qaIv) oc 1/n^ and 
Ha^Er) = E\{Er), see eq. (2.5). 

From the measured recoil spectrum in a DD experiment, TZ{Er), we can then 
extract the DM velocity distribution by 


m^ruA d ( TI{Er) \ 
Pxv‘^gA{v)dv \hA{ER)) 


(B.2) 
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Using this in the expression for the capture rate in the Sun (eqs. (3.1) and (3.2)) we 
find that the bound eq. (4.3) becomes 

_^ r^Sun r^cross 

C*Sun > 47r / drr'^pAir) / dv 

(B.3) 
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where Add indicates the mass number of the nucleus in the DD experiment, whereas 
the sum over A runs over the elements in the Sun, + 'Uesc(^)) 

rEina,^{v) 

HA{ER,r)= hA{ER)dER. (B.4) 

d EniiniA 

Eq. (B.3) corresponds to the lower bound on the capture rate if the scattering cross 
section can be factorized according to eq. (B.l). For certain models such a factorization 
may not be possible. Generalizing our bound to those cases is beyond the scope of this 
paper and we leave it for future work. 
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